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INTRODUCTION 

The  oil  embargo  of  the  1970's  revealed  the  dependence  of  the 
industrial  world  on  a  potentially  unstable  source  of  petroleum 
fuels.  Since  that  time,  increased  emphasis  has  been  placed  on 
alternate  and  renewable  fuel  sources  to  lessen  the  impact  of 
another  disruption  of  petroleum  supplies. 

One  of  the  sources  of  these  alternate  fuels  that  is 
currently  being  studied  and  developed  includes  biogas  produced  by 
anaerobic  digestion  of  livestock  wastes  and  crop  residues. 

Numerous  papers  have  been  written  on  the  theory  and  aotual 
operation  of  anaerobic  digestion  and  anaerobic  digesters  which 
use  mesophilic  and  thermophilic  bacteria  to  break  down  the  com- 
plex wastes  and  residues  to  form  methane  and  carbon  dioxide. 
Comparatively  few  studies  available  have  researched  the  use  of 
the  digester  gas  as  an  alternate  fuel  for  internal  combustion 
engines . 

Anaerobic  digestion  of  agricultural  wastes  and  residues  can 
provide  relatively  high  energy  gas  for  heating  and  fueling  pur- 
poses. Engines,  utilizing  digester  gas,  can  be  used  to  power 
other  systems  such  as  irrigation  pump  engines  or  engine  genera- 
tors to  provide  electrical  power  for  a  variety  of  uses. 

Previous  studies  have  shown  that  while  mobile  vehicles  could 
use  the  biogas  as  a  fuel,  most  vehicles  studied  laoked  adequate 
storage  capacity  for  extended  operation.   For  this  reason   it   is 


assumed   that   most   applications   using  the  digester  gas  will  be 
stationary  units. 


The  purpose  of  this  study  was  to  determine  the  necessary 
engine  modifications  to  optimize  engine  performance  and  to  iden- 
tify and  solve  some  of  the  problems  associated  with  using  the 
digester  gas  as  an  engine  fuel. 


REVIEW  OF  LITERATURE 

Methane  as  an  Engine  Fuel 

Neyeloff  and  Gunkel(1975)  reported  that  the  use  of  methane 
as  a  fuel  for  Internal  combustion  engines  has  several  advantages. 
It  has  a  high  octane  rating  which  gives  excellent  anti-knock 
qualities  as  well  as  providing  for  more  efficient  combustion  in 
engines  if  the  compression  ratio  is  increased.  Gaseous  fuels 
generally  cause  minimum  engine  carbon  deposits.  The  gaseous 
fuels  also  mix  more  thoroughly  with  air  to  burn  more  completely 
than  liquid  fuels.  Other  advantages  include  a  small  amount  of 
contaminating  pollutants,  less  sludge  in  the  lubricating  oil,  no 
wash  down  of  cylinder  wall  lubrication  during  engine  starting,  no 
tetra-ethyl  lead  to  foul  spark  plugs  and  other  engine  parts,  a 
nearly  homogeneous  mixture  in  the  cylinder,  and  less  valve  burn- 
ing. 


Flammability  limits  of  methane  in  air  are  5.3-11**  by  volume 
with  the  ideal  stoichiometric  mixture  of  about  9 . 1%  •  Neyeloff  et 
al.  (1975)  reports  that  the  ignition  limits  of  CH^  in  air  are 
affected  by  carbon  dioxide  dilutions  and  that  under  atmospheric 
conditions  the  mixture  will  not  combust  if  the  amount  of  C02  is 
greater  than  three  times  the  amount  of  methane.  The  ignition 
limits  should  be  affected  very  little  by  the  temperature  and 
pressures  developed  in  an  engine.  Yet  it  was  found  experimen- 
tally that  ignition  occurred  at  fuel  to  air  ratios  of  0.065  to 
0.185,  both  limits  higher  than  those  reported  at  atmospheric  con- 


ditions. 

Neyeloff  et  al .  (1975)  reported  that  an  ignition  advance  of 
30  degrees  before  top  dead  center  (BTDC)  gave  the  highest  effi- 
ciency for  their  CFR  engine  operated  at  900  rpm.  They  did  state 
that  at  various  times  during  the  test  runs  that  ignition  timing 
was  slightly  varied  and  little,  if  any,  changes  in  output  were 
found.  They  reported  that  specific  power  output  peaked  at  a  15:1 
compression  ratio  but  that  at  that  compression  ratio  and  above, 
audible  knocking  occurred.  The  specific  power  output  was  the 
highest  at  a  fuel  to  air  ratio  of  0.10.  A  comparison  was  made 
using  two  compression  ratios  and  different  methane  dilutions.  If 
a  compression  ratio  of  15:1  and  percent  dilution  of  0%  were  given 
a  score  of  100  for  specific  power  output,  then  a  15:1  compression 
ratio  using  a  50*  dilution  (67*  methane  and  33*  carbon  dioxide) 
would  have  a  score  of  79.  Compression  ratios  of  7.5:1  with  0* 
and  50*  dilutions  would  get  scores  of  73  and  60  respectively  for 
specific  power  output. 


Clark,  Koelsch,  Whitmar,  and  Walawender  (1978)  conducted  an 
economic  analysis  on  the  use  of  digesters  to  power  irrigation 
engines.  Advantages  of  the  digester  system  were  almost  100* 
recovery  of  residue  nutrients,  production  of  relatively  high 
energy  gas,  and  negligible  environmental  impact.  Disadvantages 
inoluded  1.)  a  need  to  distribute  about  90*  of  the  original 
solids  remaining  after  treatment,  2.)  digesters  are  not  suitable 
for  fluctuating  load  demands  without  storage,  and  3.)  digesters 
are  not  suitable  for  rapid  starts  or   shut-downs.    Additionally, 


digestion  is  subject  to  a  variety  of  instabilities  as  a  result  of 
temperature  change,  pH  variations,  carbon-nitrogen  ratio,  feed 
rate  changes,  and  the  presence  of  toxic  materials. 

The  study  was  based  on  the  operation  of  a  digester  for  a 
three-month  irrigation  season.  The  economics  improved  slightly 
if  the  digester  was  utilized  in  crop  drying  in  addition  to  irri- 
gation. Clark  et  al .  (1978)  ooncluded  that  due  to  the  high  cost 
of  gas  storage  facilities,  a  larger  digester  would  be  less  expen- 
sive and  that  any  excess  gas  produced  could  be  flared. 

Tests  conducted  by  Stahl,  Fischer,  and  Harris  (1982)  used  a 
heat  recovery  system  on  an  engine  to  heat  water  for  the  digester 
and  for  alcohol  production.  It  was  found  that  the  thermal  output 
of  the  engine  is  strongly  dependent  upon  engine  power  output. 

Engine  Performance 

Low  Compression  Spark  Ignition 

Stahl,  Harris  and  Fischer  (1982)  operating  a  3.6  liter 
engine  at  1260  rpm  produced  21  kW  using  digester  gas  with  a  55% 
methane  content.  Minimum  brake  specific  fuel  consumption 
occurred  at  a  manifold  vacuum  of  5  cm  Hg  and  an  equivalence  ratio 
in  the  range  of  0.6  to  0.8  while  peak  power  was  attained  at  an 
equivalence  ratio  of  0.825  and  a  manifold  vacuum  of  5  cm  Hg.  As 
the  throttle  was  closed,  maximum  power  for  the  throttle  setting 
was  achieved  at  equivalence  ratios  approaching  1.0.  As  the 
throttle  was  closed,  the  spark  was  advanced   for   optimum   engine 


performance.  The  researchers  also  found  that  as  the  fuel-air 
mixture  became  richer,  the  spark  must  be  retarded  from  a  maximum 
of  50°  before  top  dead  center(BTDC)  to  H5°  or  even  H0°   BTDC. 

In  two  direct  comparison  tests  performed  by  Stahl,  Harris, 
and  Fischer  (1982),  the  engine  had  slightly  higher  efficiency 
when  burning  biogas  coming  from  compressed  storage  tanks  than 
biogas  directly  from  the  digester  due  to  the  lower  water  content 
of  the  compressed  gas. 

Koelsch  (1982)  reported  that  for  an  engine  supplying  power 
to  a  generator  with  an  equivalence  ratio  greater  than  1.2  that 
misfiring  became  a  problem  while  with  an  equivalence  ratio  below 
0.78  that  the  desired  electrical  output  could  not  be  maintained. 
Electrical  efficiency  peaked  between  equivalence  ratios  of  O.78 
to  1.0.  An  equivalence  ratio  of  approximately  0.85  provided  the 
most  efficient  generation  of  electricity. 

The  test  engine  used  by  Koelsch  et  al .  (1982)  was  run  for 
1220  hours  on  unscrubbed  biogas  with  an  average  H  S  content  of 
0.4*.  The  high  sulfur  content  caused  rapid  deterioration  of  the 
buffering  capacity  of  the  oil  resulting  in  short  oil  change 
intervals.  The  capacity  of  the  oil  to  neutralize  the  effects  of 
sulfur  in  fuel  is  indicated  by  the  Total  Base  Number  (TBN).  Ini- 
tially an  oil  with  a  TBN  of  6  was  selected,  but  after  9*»  hours  of 
operation  the  TBN  was  0.46,  which  was  far  below  the  condemnation 
level  of  2.0.  The  use  of  oils  with  a  TBN  of  10  was  not  suffi- 
cient  for   unscrubbed   biogas,  slipping  below  2.0  at  55  hours  of 


operation  and  0.46  TBN  at  208  hours. 

After  1220  hours  of  operation,  all  four  rod  bearing  sets 
displayed  pitting  while  one  of  the  bearings  exhibited  severe  pit- 
ting of  the  bearing  surface  and  flaking  out  of  the  bearing  sur- 
face material.  The  pitting  and  flaking  was  attributed  to  corro- 
sive attack  on  the  bearing  surfaces  by  excessive  acidity  in  the 
oil.  An  oil  change  interval  of  250  hours  was  used  once,  while 
the  other  oil  changes  were  done  at  150  hours  or  less. 

Another  problem  was  a  slightly  excessive  wear  rate  in  some 
of  the  valve  guides,  which  was  attributed  to  the  lack  of  lubrica- 
tion in  the  fuel  and  the  lack  of  oil  reaching  the  area  between 
the  valves  and  guides.  The  other  engine  components  were  reported 
to  be  in  excellent  condition  and  were  reinstalled  in  the  engine. 


The  original  spark  plugs,  with  a  nickel  alloy  electrode, 
decayed  rapidly,  increasing  the  plug  gap  from  0.030  in.  to  0.045 
in.  in  100  hours  of  operation.  Plugs  with  an  inconel  electrode 
base  were  installed  and  after  400  hours  of  engine  operation,  the 
gap  had  inoreased  a  maximum  of  0.005  in.  Five  hundred  hour  plug 
change  intervals  were  used  with  no  apparent  problems.  Spark  plug 
selection  had  a  minor  effect  on  engine  operation  with  hotter 
plugs  to  be  somewhat  preferred.  Standard  plug  gaps  of  0.020  to 
0.030  inch  give  satisfactory  performance. 

All  spark  ignition  engines  exhibit  dramatic  reductions  in 
fuel  eoonomy  under  rich  fuel-air  mixtures  or  part  load  condi- 
tions.  At  full  load,  an  additional  0.21  m3  (7.5  ft3)   of   biogas 
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was   needed   to   produce   one   kilowatt  hour  of  electricity  at  an 
equivalence  ratio  of  1.1  as  compared  to   0.85.    Koelsch   et   al. 

•3 

(1982)  reported  that  their  unit  consumes  between  0 . 7 ^  and  0.77  oJ 
(26  and  27  ft3)  of  biogas  per  kW-hr  of  electricity  produced. 

Recommendations  for  lubrication  oils  for  digester  gas 
engines  include  oils  with  a  high  TBN  rating.  There  are,  however, 
some  trade-offs  to  consider  since  high  TBN  oils  also  have  high 
ash  contents  which  can  cause  another  set  of  engine  problems.  It 
was  also  recommended  that  engine  oil  and  coolant  temperatures  be 
maintained  above  210°  F  to  prevent  moisture  condensation  in  the 
oil  or  on  the  cylinder  walls. 

According  to  Koelsch  et  al .  (1982)  spark  timing  for  an  1800 
rpm  generator  engine  should  be  40°  BTDC  or  greater,  while  at  part 
loads  spark  advance  may  be  retarded.  It  should  not,  however,  be 
retarded  beyond  30°  BTDC.  Koelsch  (1981)  also  indicated  that  the 
spark  advance  values  given  were  the  minimum  advances  to  achieve 
maximum  power.  It  was  indicated  that  advances  of  perhaps  an  addi- 
tional 15°  showed  little  change  in  power  output.  These  values  do 
seem  to  be  the  reverse  of  the  manner  that  most  engines  with 
vacuum  advanoe  are  set  up,  indicating  that  special  distributors 
may  be  needed  for  digester  gas  engines.  A  vacuum  advance  was  not 
used  in  this  test. 


High  Compression  Spark  Ignition 

Persson  and  Bartlett  (1981)  state  that  Penn  State  had   built 
a   high-compression  spark  ignition  engine  for  biogas  by  replacing 


the  injectors  of  a  diesel  engine  with  spark  plugs  and  modifying 
the  pistons.  Preliminary  tests  on  natural  gas  indicated  higher 
thermal  efficiency  for  this  engine  than  for  lower  compression  and 
dual-fueled  engines.  The  efficiency  was  close  to  that  of  the 
diesel  fueled  engine  while  the  maximum  available  power  increased 
considerably  over  the  diesel  version. 


Dual  Fuel  Compression  Ignition 

Persson  et  al .  (1981)  report  that  diesel  engines  can  be  used 
for  biogas  with  modifications  if  some  diesel  fuel  is  used  for 
ignition  purposes.  The  original  injection  system  can  be  main- 
tained as  well  as  the  diesel's  high  compression  ratio  if  the 
engine  is  dual  fueled. 

A  constant  minimum  amount  of  liquid  fuel  must  be  injected 
for  ignition,  independent  of  the  power  output  of  the  engine.  If 
the  injection  timing  was  advanced  from  its  normal  17°  advance  to 
2k°  advance  an  increase  in  efficiency  was  obtained.  The  thermal 
efficiency  for  the  dual-fuel  engine  was  less  than  for  diesel 
fuel. 

The  normal  minimum  amount  of  diesel  fuel  that  could  be  used 
before  irregular  ignition  occurred  was  constant  for  all  tests  and 
independent  of  load  and  injection  timing.  The  amount  of  diesel 
fuel  needed  for  proper  ignition  depends  on  the  design  of  the  fuel 
system  and  the  engine  design.  At  100$  load  and  maximum  effi- 
ciency, the  liquid  fuel  delivered  less  than  10$  of  the  necessary 
fuel  energy  while  at  10$  load,  the  diesel   supplied   40$   of   the 
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energy.  Since  it  is  desirable  to  use  mainly  the  biogas  for 
power,  the  engine  should  operate  as  close  to  full  load  as  possi- 
ble. Persson  et  al.(198l)  also  reported  that  the  dual-fuel 
engine  operated  quieter  than  did  the  100*  diesel  engine. 

Persson  and  Bartlett  (1981)  report  that  some  of  the  water 
vapor  in  the  gas  should  be  removed  to  prevent  condensation  prior 
to  use  as  fuel,  as  this  will  aid  in  preventing  fouling  of  the  gas 
lines  and  valves.  It  may  also  be  desirable  to  maintain  a  fairly 
high  engine  temperature,  even  when  not  in  operation  to  prevent 
water  vapor  condensation.  Reductions  of  H_S  concentrations  to 
approximately  1  mg  H  S  per  liter  was  also  recommended.  It  was 
noted  that  high-strength  steel  seems  to  be  more  susceptible  to 
hydrogen  sulfide  attack  than  cast  iron  or  lower  grade  steels. 

Ortiz-Canavate,  Hills,  and  Chancellor  (1981)  used  a  conven- 
tional diesel  engine  as  a  dual-fuel  engine  using  diesel  and  syn- 
thetic biogas  (natural  gas  and  carbon  dioxide  mixtures). 
Twenty-three  degree  injector  advance  instead  of  the  the  normal 
19°  advance  seemed  more  favorable  to  avoid  knocking  and  misfiring 
at  low  torque  levels  or  high  engine  speeds. 

The  air  ratio,  the  actual  air  used  to  stoichiometric  ratio, 
was  1.2:1  for  dual  fuel  as  compared  to  1.5:1  for  diesel.  Effi- 
ciency of  the  dual  fuel  engine  when  operating  at  high  torque  lev- 
els and  low  engine  speeds  was  equivalent  to  that  of  the  engine 
operating  solely  on  diesel  fuel.  At  higher  speeds,  the  air-fuel 
mixture  was  less  than  stoichiometric  and  efficiency  dropped  while 
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at  low  torques,  the  efficiency  of  the  dual-  fuel  was   much   lower 
than  the  diesel  and  the  engine  did  not  run  smoothly. 

Digester  Gas  Clean-up 


Digester  gas  cleanup  is  important  from  several  aspects.  It 
may  be  necessary  to  treat  the  gas  to  meet  pipeline  specifications 
if  the  gas  is  to  be  used  commercially.  Pipeline  specifications 
are  important  from  the  standpoints  of  uniform  heating  value,  cor- 
rosiveness,  and  efficiency  of  transportation  and  heating.  Diges- 
ter gas  usually  will  have  four  constituents;  methane,  carbon 
dioxide,  hydrogen  sulfide,  and  water  vapor.  Of  these  the  hydro- 
gen sulfide  and  the  water  vapor  cause  serious  corrosion  problems 
for  pipelines,  valves,  furnaces,  and  engines.  These  constituents 
should  be  removed  as  soon  as  is  feasible.  The  carbon  dioxide  is 
generally  non-corrosive  except  when  absorbed  in  water. 

Dynatech  (1978)  conducted  a  study  of  different  systems  for 
the  removal  of  water  vapor,  HO  and  C0?.  The  study  was  conducted 
for  three  raw  gas  flowrates  and  for  two  end  delivery  pressures. 
The  flowrates  were  100,000  standard  cubic  feet  per  day  (SCFD),  1 
million  (MM)  SCFD,  and  3MMSCFD;  the  0.1MMSCFD  corresponds  to  the 
expected  output  from  a  2000  head  beef  feedlot.  The  end  pressure 
of  the  sweet  (clean)  gas  was  125  psi  for  intrastate  pipeline 
transmission  and  1000  psi  for  interstate  pipelines.  No  credit 
was  assumed  for  sale  of  carbon  dioxide  or  hydrogen  sulfide, 
although  it  was  estimated  that  as  muoh  as  70  tons  per  day  of  car- 
bon dioxide  could   be   recovered   from   the   3MMSCFD   flow.    The 
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recovery  of  the  hydrogen  sulfide  was  felt  to  be  economically 
unfeasible  with  production  of  4.5,  45,  and  134-lbs.  per  day  from 
the  three  flow  rates  respectively. 

Pipeline  specifications  require  that  the  gas  be  at  least  97* 
methane,  less  than  3  mole  percent  carbon  dioxide,  less  than  0.25 
grains  per  100  SCF  (0.00036$)  H  s,  and  water  vapor  shall  not 
exceed  seven  pounds  per  MMSCF.  The  heating  value  of  the  gas 
shall  not  be  less  than  950  BTU  per  SCF. 

After  examining  the  expected  costs  of  different  systems  now 
available,  Dynatech  decided  that  the  water  scrubbing  process  for 
carbon  dioxide  removal,  hydrogen  sulfide  removal  with  a  process 
developed  by  Eickmeyer  &  Associates,  and  water  vapor  removal  with 
triethylene  glycol  (TEG)  would  be  the  most  economical.  Many  of 
the  systems  now  used  by  the  commercial  gas  industry  were 
developed  for  treatment  of  much  larger  quantities  of  gas  and  were 
therefore  not  economical  to  be  scaled  down  for  use  of  biogas 
cleaning. 

For  the  Dynatech  system,  the  costs  for  the  three  flow  rates, 

0.1,   1,   and  3MMSCFD,  were  $1,943,  $1,169,  and  $.780  per  million 

BTU  at  130  psl  delivery  pressure  and  $1,982,  $1,178,  and  $.790 
respectively,  per  million  BTU  at  1020  psi  delivery. 

The  water  scrubbing  process  has  been  used  for  centuries,  and 
uses  water  to  absorb  carbon  dioxide.  The  process  works  best  at 
high  pressures  (above  300  psia  is  recommended)  and  at  low  tem- 
peratures.   It  is  quite  dependent  on  the  partial  pressure  of  the 
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CO-.  Dynatech  suggested  the  use  of  a  lagoon  to  desorb  the  carbon 
dioxide  at  atmospheric  pressures  while  other  schemes  have  used  a 
series  of  flash  desorbers  so  that  the  C0?  may  be  recovered.  If 
carbon  dioxide  recovery  is  desired,  there  are  other  commercially 
available  methods  that  are  more  efficient. 

The  Eickmeyer  process  for  hydrogen  sulfide  removal  uses  an 
undisclosed  proprietary  chemical  solvent.  The  sour  gas  is  bub- 
bled through  the  solvent  in  a  shallow  tank  at  atmospheric  pres- 
sure with  almost  complete  removal  of  H  S.  The  solvent  is  highly 
selective  for  H  s  The  amount  of  solution  is  on  the  order  of  1000 
gal/MMSCFD  of  gas  flow.   The  process  is  non-regenerative. 


Water  absorption  uses  either  diethylene  glycol  (DEG)  or 
triethylene  glycol  (TEG).  Advantages  of  this  process  are  unusual 
hygroscopicity  of  the  solution,  good  stability  to  heat  and  chemi- 
cal decomposition,  low  vapor  pressures,  and  ready  availability  at 
moderate  cost.  In  a  typical  unit,  water  vapor  is  continuously 
absorbed  from  the  process  gas  stream  by  countercurrent  contact 
with  a  highly  concentrated  glycol  stream  (95  to  99*)  in  a  packed 
or  bubbly  tray  column.  Dynatech  reported  that  that  they  based 
their  review  on  a  system  that  had  a  TEG  concentration  of  99.8 
weight  percent  for  low  pressures  and  for  high  pressure,  a  TEG 
concentration  of  98.7  weight  percent.  Regeneration  of  the  rich 
glycol  solution  is  usually  accomplished  by  inert  gas  stripping 
with  the  application  of  heat.  When  maximum  dehydration  is 
required  in  large  Installations,  vacuum  regeneration  can  also  be 
used . 
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INVESTIGATION 
Objectives 
The  objectives  of  the  investigation  were: 

1  .   To  compare  the  performance  of  the  engine  operation  on 
two  different  sources  of  anaerobic  digester  gas  with 
natural  gas. 

2.  To  determine  the  effect  of  the  fuels  on  the  engine 
components . 

3.  To  determine  the  changes  in  spark  advance  needed  to 
optimize  operation. 

Theory 

Patterson  et  al .  (1972)  report  that  the  combustion  process 
occurs  in  two  steps,  the  ignition  delay  and  the  pressure-rise 
period.  The  ignition  delay  period  is  the  time  required  for  a 
fuel  and  air  mixture  to  auto-ignite  when  subjected  to  a  high  tem- 
perature at  a  given  pressure.  The  length  of  the  ignition  delay 
period  depends  upon  charge  pressure,  fuel-air  ratio,  and  octane 
number  of  the  fuel.  Ignition  delay  is  independent  of  engine 
speed;  thus,  at  high  engine  speeds  the  ignition  delay  occupies  a 
larger  number  of  crank  angle  degrees  than  at  lower  speeds. 
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The  pressure-rise  period  depends  on  the  type  of  fuel  and  the 
turbulence,  among  other  factors.  Since  the  flame  velocity  is 
proportional  to  the  turbulence  which  in  turn  is  proportional  to 
engine  speed,  the  pressure-rise  period  occupies  about  the  same 
number  of  crank  angle  degrees  at  different  speeds. 

Methane  is  a  slow  burning  fuel,  therefore,  the  timing 
advance  of  engines  fueled  with  methane  needs  to  be  increased  to 
maximize  the  engine  output  power.  If  the  spark  advance  is  not 
advanced  enough,  the  fuel  will  not  have  burned  completely  before 
being  exhausted  through  the  exhaust  valves.  This  exposes  the 
exhaust  valves  and  manifold  to  higher  temperatures  and  thus  the 
expected  life  of  these  components  is  shortened.  If  the  timing  is 
too  advanced,  the  peak  pressure  rise  will  occur  before  top  dead 
center.  Knocking  may  occur  and  performance  suffer  as  a  result. 
Dilution  of  methane  with  inert  gasses  such  as  carbon  dioxide 
reduces  the  flame  velooity  of  the  methane  mixtures.  This  dilu- 
tion also  causes  the  need  for  greater  spark  advance. 

Materials  and  Experimental  Equipment 

Test  Fuels 

The  fuels  tested  were  natural  gas,  anaerobic  digester  gas 
derived  from  beef  manure,  and  anaerobic  digester  gas  derived  from 
swine  manure.  The  natural  gas  was  supplied  by  Kansas  Power  and 
Light  Company  at  Kansas  State  University  in  the  Agricultural 
Engineering's  Engine  Test  Laboratory.  The  beef-based  biogas  was 
produced   at   the   D.S.D.A.   Roman  L.  Hruska  Meat  Animal  Research 
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Center  at  Clay  Center,  Nebraska  while   the   swine-based   gas   was 
made  available  at  the  Del  Valle  Hog  Farm  at  Austin,  Texas. 

The  biogas  at  both  locations  was  characterized  by  a  50  to 
65$  methane  content,  with  carbon  dioxide  concentrations  of  35  to 
50$;  the  remainder  of  the  constituents  being  hydrogen  sulfide  and 
water  vapor. 

Test  Equipment 

Test  Engine  The  test  engine  was  a  Ford  industrial  engine,  Model 
KSG-411.  The  engine  was  equipped  with  an  Impoo  CA50-510  carbure- 
tor designed  for  operation  on  natural  and  low  energy  gas.  The 
engine  was  originally  equipped  with  a  larger  carburetor  to  allow 
for  a  greater  amount  of  low  energy  biogas  to  enter  the  engine. 
Engine  performance  was  not  up  to  expectations  and  analysis  of  the 
exhaust  gasses  revealed  high  oxygen  and  low  carbon  monoxide  lev- 
els indicating  a  very  lean  fuel-air  mixture.  The  original  car- 
buretor was  not  properly  matched  to  the  engine.  It  was  designed 
for  a  larger  displacement  engine.  The  carburetor  fuel  Jets  had 
the  capability  of  providing  adequate  gas  for  the  small  engine  but 
a  large  air  flow  was  needed  to  provide  the  pressure  drop  across 
the  air  venturi  to  obtain  the  respective  gas  flow.  This  resulted 
in  poor  performance  from  a  lean  mixture.  A  solution  to  this 
problem  proved  to  be  a  smaller  carburetor  with  increased  gas  sup- 
ply pressure  to  supply  the  needed  quantity  of  hydrocarbons  per 
given  volume.   The  engine  specifications  are  shown  in  Table  1. 
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n 

- 

80.98 

mm 

53.29 

mm 

1  .1 

liters 

8:1 

- 

1500-: 

3600 

rpm 

18 

kW 

Table  1.  Ford  KSG-411  Test  Engine  Specifications 

Specification Value      Dnits 

Cylinders 

Bore 

Stroke 

Displacement 

Compression  Ratio 

Rated  Speed 

Power ( continuous  on  gasoline) 


The  engine  is  equipped  with  a  fan  and  radiator,  alternator, 
and  a  muffler.  A  Rockwell  PTO  clutch  and  an  instrument  panel  were 
installed.  The  panel  featured  a  starter  key  switch,  a  coolant 
temperature  gauge,  an  oil  pressure  gauge,  a  voltmeter,  and  a  ver- 
nier type  throttle  control. 

Dynamometer  ajil  Instrumentation.  A  model  70523  Cessna  pump 
equipped  with  a  Parker-Hannif in  R10PH-11-BL  pressure  control 
valve  was  used  as  a  dynamometer  load  for  all  the  tests.  A  90.7 
kg  (200-lb.)  Lebow  3167  load  cell  installed  on  an  279-4  mm  lever 
arm  was  used  to  determine  torque.  The  load  cell  signal  was 
received  by  a  Daytronic  3270  strain  gage  conditioner  and  indica- 
tor. Engine  speed  was  determined  by  a  60-tooth  gear  and  an  Elec- 
tro 3010AN  magnetic  pickup,  with  the  signal  received  by  a  Day- 
tronic 3240  frequency  conditioner  and  indicator. 

The  digester  gas  was  filtered  through  two  Winslow  gas 
filters  (Model  981136-B)  arranged  in  parallel  to  avoid  overload- 
ing the  filters.  A  Rockwell  415  gas  meter  and  indicator  was  used 
in  the  first  series  of  tests  to  provide  a  visual  indication  of 
gas  consumption.   A  magnetic  reed  switoh  was  used   to   provide   a 
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pulse  signal  to  the  data  acquisition  system  to  determine  the  rate 
of  gas  use. 

Copper-constantan  thermocouples  were  used  to  measure  the 
temperatures  of  the  oil,  the  fuel,  the  inlet  and  outlet  coolant, 
and  ambient  intake  air.  The  exhaust  gas  temperature  was  measured 
with  an  iron-constantan  thermocouple. 

For  the  second  series  of  tests,  absolute  fuel  and  atmos- 
pheric air  pressures  were  measured  by  Setra  Systems  Inc.  Model 
204  0-315  kPa  (0-50  psi)  variable  capacitance  absolute  pressure 
transducers . 

Air  flow  measurement  used  a  31.75  mm  (1.25  in.)  diameter 
ASME  nozzle  mounted  in  a  114  liter  (30  gal.)  plenum  chamber.  The 
large  plenum  chamber  was  needed  to  smooth  the  pulses  caused  by 
the  cyclic  intake  of  air.  The  pressure  drop  across  the  nozzle 
was  measured  with  a  Setra  Systems  Ino.  Model  239  0-1.4  kPa  (0-0.2 
psi)  variable  capacitance  differential-pressure  transducer.  The 
fuel  flow  measurement  Involved  measuring  the  pressure  drop  across 
a  12.7  mm  (0.5  in.)  diameter  ASME  nozzle.  This  pressure  drop  was 
also  measured  with  a  Setra  Systems  Inc.  Model  239  differential- 
pressure  transducer.  The  differential  transducers  were  cali- 
brated against  a  Meriam  Model  34FB2  micromanometer .  The 
discharge  coefficient  for  the  flow  nozzles  was  determined  by 
using  a  relationship  proposed  by  Benedict  (1966),  relating  the 
most  probable  disoharge  coefficient  to  the  Reynold's  number.  The 
fuel  and  air  supply  system  is  illustrated  in  Figure  1 . 
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Figure  1.   Fuel  and  Air  Supply  Instumentation 
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The  gas  and  air  measurement  system  used  performed  well  at 
low  speeds,  but  at  engine  speeds  of  approximately  2000  rpm  a 
resonance  frequency  apparently  appears  and  the  indicated  gas  and 
air  use  is  much  lower  than  is  realistic.  Specific  fuel  consump- 
tion and  thermal  efficiencies  were  unreasonably  low  and  high, 
respectively,  and  will  not  be  included  in  the  results. 

The  data  acquisition  system  used  a  Synertek  Systems  Corp. 
SYM  1  computer  to  collect  data.  The  test  data  was  stored  on 
cassette  tape  until  the  data  could  be  processed  on  a  DEC  PDP  11- 
31  computer. 

The  SYM  1  featured  a  16-channel,  analog-to-digital  board 
with  two  pulse  counting  channels  to  facilitate  data  collection. 
Data  was  taken  at  15-second  intervals  and  stored  in  three  minute 
blocks  which  could  be  averaged  using  Fortran  software.  The  For- 
tran program  is  inoluded  in  Appendix  B. 

Test  Procedures 

An  engine  break-in  period  of  50  hours  was  used  prior  to 
operating  the  engine  on  biogas.  A  90-minute  break-in  cycle  pro- 
posed by  the  Engine  Manufacturer's  Association  (EMA)  to  the 
United  States  Department  of  Agriculture  was  used  with  commercial 
natural  gas  to  fuel  the  engine.  The  break-in  cycle  was  as  fol- 
lows : 

1.  10  minutes  at  low  idle 
(1200  rpm) 

2.  10  minutes  at  1/2  rated  speed  and  no  power 
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(  1440  rpm) 


15  minutes  at  3M  rated  speed  and  1/2  rated  power 
(2160  rpm  and  17.6  N-m) 


55  minutes  at  rated  speed  and  rated  power 
(2880  rpm  and  35  N-m) 


The  actual  engine  tests  consisted  of  a  200-hour   cycle   also 

recommended  by  EMA .   The  test  cycle  was  as  follows: 

1  .  30  minutes  at  low  idle 
(  1 1400  rpm) 


2.  60  minutes  at  rated  speed  and  rated  power 
(2880  rpm  and  35  N-m) 


3.  60  minutes  at  peak  torque 
(2400  rpm  and  36.9  N-m) 


4.  30  minutes  at  80$  rated  speed  and  25%    rated  power 
(2300  rpm  and  11.6  N-m) 


For  the  first  series  of  tests,  due  to  the  availability  of 
the  biogas,  two  three-hour  cycles  were  run  per  day  with  an  aver- 
age of  25  hours  of  engine  operation  per  week.  Torque  curves  were 
run  every  25  hours  and  performance  was  mapped  at  the  beginning, 
middle  and  end  of  the  200  hour  test.  Data  was  taken  at  200  rpm 
intervals  and  values  for  the  performance  map  at  intervals  of 
one-tenth  full  torque.  The  engine  was  run  twenty-five  hours  on 
natural  gas  prior  to  the  second  test  to  purge  the  engine  of  the 
first  test's  residuals.  During  the  second  test  series,  the 
engine  was  run  for  five  three-hour  loading  cycles  for  a  total  of 
fifteen  hours  per  day.   The  two-hundred  hour  test   was   completed 
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in  approximately  two  weeks.   Collection  of  data  was  conducted  the 
same  as  for  the  first  series  of  tests. 

Oil  samples  were  collected  at  25-hour  intervals  and  sent  to 
Farmland  Industries  to  check  for  wear  metal  levels  and  oil 
deterioration.  The  oil  and  filters  were  changed  at  100-hour 
intervals . 

flesvltg 

Engine  Performance  on  Natural  Gas 

The  engine  power  peaked  at  16.0  bkW,  89.5%  of  the  rated 
power  of  the  engine  on  gasoline.  Peak  power  occurred  at  3600 
rpm.  The  peak  torque,  54.2  N-m,  occurred  at  2400  rpm.  The 
torque  and  power  curves  are  illustrated  for  reference  to  the  bio- 
gas  tests  in  Figures  2  through  5. 

Test  Series  1:   Engine  Performance  on  Beef  Manure  Biogas 

Eight  torque  curves  were  conducted  while  operating  on  the 
beef-based  biogas.  The  curves  varied  depending  on  the  biogas 
properties  at  that  time.  The  methane  content  was  usually  about 
53  percent  although  fluctuations  of  up  to  7  percent  were  seen  in 
successive  weeks.  The  gas  pressure  also  was  subject  to  small 
variations.  Average  torque  and  power  curves  were  calculated  from 
the  eight  torque  curves.  The  results  are  plotted  in  Figures  2 
and  4,  respectively,  as  are  the  95f  confidence  limits  for  the 
expected  limits  of  day  to  day  operation.  The  peak  power  of  the 
average   curve   was   12.3  bkW,  77%    of  the  peak  natural  gas  power. 
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However,  the  highest  observed  power  was  13-3  bkW  }  83$  of  natural 
gas  power.  The  highest  value  on  the  average  torque  curve  was 
38.7  N-m  which  occurred  at  2000  rpm  while  the  highest  observed 
torque,  11. 45  N-m,  was  at  2200  rpm. 

Test  Series  2: Engine  Performance  on  Swine  Manure  Biogas 

Eight  torque  curves  were  run  at  25-hour  intervals  while  the 
engine  was  fueled  by  swine-derived  digester  gas.  The  methane 
content  of  the  gas  was  consistently  64  percent  with  only  occa- 
sional slight  variation  of  a  maximum  two  percent  change.  Average 
torque  and  power  curves  with  95$  confidence  limits  were  calcu- 
lated from  the  eight  torque  curves.  The  results  are  plotted  in 
Figures  3  and  5.  The  peak  power  of  the  average  curve  was  15.0 
bkW,  94%  of  the  peak  natural  gas  power  and  122$  of  the  peak  aver- 
age beef-derived  gas.  The  highest  power  observed  was  15.2  bkW, 
95$  of  the  peak  on  natural  gas.  The  highest  torque  of  the  aver- 
age curve  was  41.6  N-m  which  occurred  at  2000  rpm  while  the 
highest  observed  torque,  43.3  N-m  was  also  at  2000  rpm.  The  peak 
on  the  average  torque  curve  was  77$  of  the  peak  on  natural  gas 
and  107$  of  the  peak  of  the  average  curve  of  the  first  test 
series  . 
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Figure  2.   Torque  Curve:  Test  Series  1 
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Figure  3.   Torque  Curve:  Test  Series  2 
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Figure  3.   Power  Curve:  Test  Series  1 
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Figure  5.  Power  Curve:  Test  Series  2 


28 


Engine  Performance 

The  hydrocarbon  content  of  the  fuel,  the  amount  of  water  in 
the  gas,  and  the  fuel  pressure  appear  to  be  the  controlling  fac- 
tors in  the  engine  performance.  The  biogas  supply  pressure  was 
quite  important.  If  the  pressure  was  less  than  750  mm  H  n,  the 
performance  was  sluggish.  At  times  in  the  first  test  series,  dur- 
ing the  3  hour  engine  cycle,  water  vapor  in  the  gas  would  reduce 
the  engine  output  enough  that  the  dynamometer  load  would  kill  the 
engine.  At  those  times  the  engine  could  operate  at  1800  rpm  with 
a  maximum  torque  setting  of  approximately  10  N-m.  There  were  no 
problems  with  water  vapor  at  the  Texas  facility.  This  can  prob- 
ably be  attributed  to  the  difference  in  the  storage  tank  size 
between  the  two  facilities.  The  Nebraska  facility  used  two  1000 
gallon  propane  tanks  to  hold  the  gas  production  while  the  Texas 
facility  used  a  60,000  gallon  butane  tank  for  storage.  The 
larger  tank  would  have  less  turbulence  allowing  the  water  vapor 
to  condense  and  settle  out  of  the  gas.  The  tank  could  then  be 
drained  of  the  condensed  water  periodically.  The  methane  content 
of  the  gas  changed  engine  performance  very  little.  A  5%  drop  of 
methane  content  gave  only  slightly  noticeable  performance 
deterioration . 

Spark  Advance 


Shown  in  Figure  6  is  the  spark  advance  of  the  engine  operat- 
ing on  digester  gas  at  full  torque.   The  spark  advance  shown  is 
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Figure  6.   Spark  Advance:  Test  Series  1  and  2 
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the  minimum  advance  for  best  torque.  The  difference  in  the  two 
curves  is  due  to  the  difference  in  methane  content  of  the  two 
biogas  supplies.  The  richer  gas  from  the  swine  has  less  carbon 
dioxide  dilution,  resulting  in  less  slowing  of  the  flame  velocity 
in  the  combustion  chamber.  If  the  advance  is  increased  beyond 
these  values,  little  change  in  power  was  detected. 

Engine  Oil  Analysis 

Oil  samples  were  collected  at  25-hour  intervals  and  were 
analyzed  by  Farmland  Industries  for  wear  metal  concentrations. 
Metal  concentrations  (ppm)  that  are  detected  are  iron,  chromium, 
aluminum,  copper,  lead,  tin,  silver,  nickel,  silicon,  sodium, 
boron,  magnesium,  calcium,  barium,  phosphorus,  and  zinc.  The  oil 
viscosity  at  100  C,  the  presence  of  water  and  antifreeze  in  the 
oil,  total  solids,  and  oxidation  were  also  monitored.  Iron  con- 
centrations indicate  cylinder  wall  wear.  Chromium  indicates  ring 
wear.  Aluminum  points  to  piston  wear.  Copper,  lead,  and  tin 
indicate  bearing  wear,  while  silicon  indicates  that  dirt  or  sand 
is  being  induced  into  the  engine  usually  from  a  leaking  air 
cleaner.  Total  Base  Number  (TBN)  was  monitored  during  the  second 
test  series.  TBN  is  an  indicator  of  the  remaining  buffering  qual- 
ities of  the  oil.  Table  2  illustrates  some  of  the  wear  metal 
concentrations  (ppm)  in  the  oil  at  25  hour  intervals. 

Condemnation  limits  are  not  given  as  variations  between 
engines  are  great;  abrupt  changes  of  metal  concentrations  between 
samples  should  be  the  primary  indicator  of  problems  and  of  exces- 
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sive  wear. 


Table  2.   Ford  KSG-411  Engine  Oil  Wear  Metal  Analysis 


Test 

Series 

1  :   Bee 

f  ; 

Manure 

Deri 

ved  Fuel 

Hours 

Iron 

Chromi 

urn 

Alumi 

num 

Copper 

Lead 

Tin 

Silicon 

TBN 

24 

11 

0 

0 

7 

3 

4 

10 

- 

50 

12 

0 

0 

8 

3 

7 

9 

- 

75 

17 

0 

0 

9 

3 

8 

10 

- 

100 

16 

0 

0 

9 

2 

5 

7 

- 

125 

7 

0 

8 

3 

1 

0 

6 

- 

150 

10 

0 

0 

4 

0 

0 

5 

- 

175 

19 

0 

0 

7 

5 

3 

7 

- 

200 

11 

0 

5 

5 

0 

0 

4 

- 

Test 

Series 

2:   Swi 

ne 

Manure 

Der 

ived  Fuel 

25 

18 

0 

7 

10 

3 

0 

6 

10.1 

50 

23 

0 

8 

15 

3 

0 

5 

10.7 

75 

34 

0 

0 

26 

0 

3 

7 

10.0 

100 

34 

0 

2 

23 

0 

9 

7 

9-9 

125 

1  1 

0 

2 

7 

0 

5 

9 

10.3 

150 

14 

0 

9 

12 

2 

0 

6 

10.4 

175 

20 

0 

7 

18 

1 

0 

7 

10.2 

200 

24 

0 

0 

17 

0 

0 

7 

9-9 

The  wear  metal  concentrations  present  in  the  oil  indicate 
that  no  significant  wear  has  occurred  in  400  hours  of  engine 
operation  on  digester  gas.  The  Total  Base  Number  after  100  hours 
of  use  is  still  quite  high  suggesting  that  the  oil  change  inter- 
val could  be  significantly  extended.  A  greater  period  between 
oil  changes  will  decrease  maintenance  time  and  cost,  improving 
the  economics  of  operating  the  engine  as  a  power  source. 

Engine  Deposits  and  Enftlne  Wear 


The  head  was  removed  after  the  first  200  hour  test  series. 
The  small  amount  of  deposits  visible  on  the  head  and  on  the  upper 
cylinder  walls  was  soft  and  easily  removed.  The  engine  was  dis- 
mantled  after   the   second   test  series  and  significant  deposits 
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were  found  on  the  head,  the  pistons,  and  the  valves.  The  depo- 
sits on  the  head  were  nearly  to  the  point  of  interfering  with  the 
valve  seatings  and  looked  as  though  valve  problems  would  have 
soon  developed.  Photographs  of  the  deposits  on  the  components 
are  located  in  Appendix  A.  An  elemental  analysis  and  an  electron 
microscope  scan  were  conducted  on  the  deposits.  The  results  show 
that  carbon,  hydrogen,  nitrogen,  and  sulfur  are  present  with 
trace  amounts  of  zinc,  calcium,  iron, 

The  journal  bearings  and  other  engine  components  were  exam- 
ined after  the  second  test  series.  No  excessive  wear  was  evi- 
dent. No  pitting  or  warnings  of  impending  failure  of  engine  com- 
ponents could  be  seen,  unlike  results  found  by  researchers  at 
Cornell.  The  main  bearings  were  plastigaged  to  determine  wear. 
The  clearances  were  well  within  new  bearing  tolerances. 

Digester  Gas  Analysis 

An  attempt  to  determine  the  effectiveness  of  the  biogas 
filters  led  to  gas  chromatography  analysis  of  the  swine  biogas. 
A  50  meter  column  was  used  with  a  mass  spectrometer  to  determine 
the  compounds  present  in  the  gas.  Compounds  with  mass  number 
larger  than  30  were  being  detected  by  the  analysis  with  only  car- 
bon dioxide  and  hydrogen  sulfide  appearing.  No  mercaptans,  that 
is  compounds  with  an  S-H  radical,  appeared  in  the  analysis  con- 
trary to  literature  written  by  the  filter  manufacturer  reporting 
the  presence  of  the  mercaptans.  Gas  samples  were  taken  at  three 
points   in   the   system.   The  first  was  taken  from  the  top  of  the 
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digester  before  any  filtering  occurred.  The  second  sample  had 
passed  through  an  iron  sponge  and  had  been  compressed  into  the 
storage  reservoir.  The  third  gas  sample  had  also  passed  though 
the  iron  sponge  and  the  compressor  as  well  as  filtering  though 
the  Nelson  filters.  The  gas  at  each  sampling  point  was  approxi- 
mately 36  percent  C02.  The  results  indicate  that  no  significant 
reduction  of  the  hydrogen  sulfide  was  taking  place  in  either  of 
the  filters.  The  tests  were  set  up  to  determine  the  quantity  of 
H_S  as  a  fraction  of  the  quantity  of  COp.  The  results  are  shown 
in  Table  3. 


Table  3.   Fraction  of  H2s  Present  in  36*  C02  Biogas 


Filtering 


$Total-C02    $Total-H2S 


none 

none 

none 

iron  sponge 

iron  sponge 

iron  sponge 

iron  sponge  and  Nelson  filter 

iron  sponge  and  Nelson  filter 

iron  sponge  and  Nelson  filter 


98.69 

1.31 

98.60 

1 .40 

98.54 

1.46 

98.96 

1.04 

98.64 

1.36 

98.84 

1.16 

99.18 

0.82 

98.39 

1 .61 

98.63 

1.37 

3* 

CONCLUSIONS 

It  is  desirable  to  remove  the  hydrogen  sulfide  from  the  bio- 
gas  as  close  to  the  digester  as  possible  to  minimize  the  damage 
caused  by  corrosion  of  metal  surfaces.  The  filters  that  were 
used  to  remove  the  hydrogen  sulfide  from  the  digester  gas  proved 
to  make  insignificant,  if  any,  reduction  of  the  hydrogen  sulfide 
content . 

Larger  storage  vessels  to  allow  water  vapor  to  condense  are 
desirable  both  for  optimizing  engine  output  and  for  minimizing 
the  corrosion  of  engine  components. 

Proper  sizing  of  the  engine,  carburetor,  and  accessories  are 
essential  for  satisfactory  performance. 

Engine  wear  is  not  considered  to  be  a  problem,  at  least  for 
lubrication  oil  change  intervals  of  100  hours.  Further  tests 
should  be  conducted  to  determine  the  maximum  period  of  safe 
engine  operation  per  oil  change. 

The  head  should  be  removed  periodically  for  inspection  to 
identify  any  problems  that  might  be  corrected  by  minor  altera- 
tions, such  as  addition  of  a  filter  to  remove  oil  leaking  out  of 
the  compresser  into  the  gas. 

Optimum  spark  advance  and  fuel  gas  to  air  mixtures  is  depen- 
dent on  the  methane  content  of  the  gas.  Small  changes  of  timing 
or  mixture  affeot  peak  torque  and  power  only  to  a  slight  extent. 
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Operation  of  an  engine  on  digester  gas  should  yield  torque 
and  power  outputs  of  80  to  95  percent  of  operation  on  natural 
gas.  This  will  depend  on  the  energy  content  of  the  biogas  gas 
which  is  dependent  upon  the  type  of  animal,  the  feed  ration,  the 
operation  of  the  digester,  and  the  digester  design. 
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This  study  shows  that  with  minor  engine  adjustment,  opera- 
tion and  performance  of  a  spark  ignition  engine  with  an  8  to  1 
compression  ratio  on  digester  gas  is  satisfactory.  Adjustments 
inolude  advancement  of  the  timing,  enrichment  of  the  gas  to  air 
mixture,  and  an  increased  fuel  supply  pressure. 

The  results  of  the  study  show  that  biogas  supplies  80  to  95 
percent  of  natural  gas  power.  Lubricating  oil  did  not  exhibit 
serious  signs  of  deterioration  with  100-hour  oil  change  inter- 
vals. 

Problems  encountered  inolude  water  vapor  and  hydrogen  sul- 
fide concentrations  in  the  gas.  In  long  term  engine  tests, 
severe  corrosion  of  engine  components  may  occur.  Engine  deposits 
were  quite  heavy  after  200  hours  of  operation  on  swine  digester 
gas. 

These  results  can  be  projected  only  for  short   term   use  of 

anaerobic   digester   gas   as  an  engine  fuel.   Long  term  tests  are 

needed  to  verify  any  assumptions  derived  by  extrapolation  of  the 
results. 
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SUGGESTIONS  FOR  FURTHER  RESEARCH 

Further  studies  should  be  made  in  the  development  and  effec- 
tiveness testing  of  filters  for  removing  hydrogen  sulfide  from 
anaerobic  digester  gas.  The  filter  needs  to  be  inexpensive  and 
have  low  maintenance  requirements. 

Research  should  be  performed  to  better  utilize  the  high 
octane  content  of  the  biogas  by  testing  high  compression,  spark 
ignition  engines  to  inorease  the  efficiency  of  the  engine  cycle. 
If  these  engines  would  also  have  increased  turbulence  in  the 
combustion  chamber,  the  flame  speed  would  probably  increase. 
This  would  allow  the  spark  advance  to  be  retarded  so  that  a 
sharper  pressure  rise  could  occur  near  top  dead  center  to  optim- 
ize the  power  output  of  the  engine. 

A  lubricating  oil  with  a  high  buffering  capacity  should 
allow  greater  time  intervals  between  oil  changes.  An  alternative 
to  this  would  be  to  develop  acid-resistant  engine  components. 

For  an  industry  that  has  demands  on  time  that  are  essen- 
tially fully  extended,  the  entire  system  needs  to  operate 
automatically  with  a  high  reliability.  This  system  would  need  a 
computer  to  monitor  and  control  the  entire  process  from  loading 
the  digester  to  running  the  engine  generator  during  peak  electric 
demand  periods.  An  operator  should  have  to  do  only  periodic 
maintenance  such  as  changing  oil  in  the  engine. 
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APPENDIX  A.   PHOTOGRAPHS  OF  ENGINE  DEPOSITS 


41 


Figure  7.   Piston  Side  View  after  200  Hours  of 
Engine  Operation,  Swine  Digester  Gas 


Figure  8.   Top  of  Piston  after  200  Hours  of 

Engine  Operation,  Swine  Digester  Gas 
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Figure  9.   Face  o_ 

Engine  Operation, 


f  Head  after  200  Hours  of 

Swine  Digester  Gas 


Figure  10.   Intake  and  Exhaust  Valves  after 
200  Hours  of  Engine  Operation, 
Swine  Digester  Gas 
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APPENDII  B. FQRIflAH  PR.QGBAMS 


c   MAIN  PROGRAM  "ENGINE. F" 
inolude  'flog.h' 


c  this  program  calculates  torque,  power,  mean  effective 

c  pressure,  fuel  flow  rate,  specific  fuel  consumption, 

c  thermal  efficiency,  correction  factor  for  temperature 

c  and  pressure  for  values  on  the  torque  curve  for  a  67 

c  cubic  inoh  displacement  Ford  engine  on  a  dynamometer 

c  with  a  load  cell  on  an  11  inoh  lever  arm. 

c  c  or  cf  prefix  indicates  corrected  for  temp,  and  press, 

c  m  suffix  indicates  metric  units 

real  pgas,  paatm,  frac 
real  tmps( 8) ,gas , load , rpm , pt ( 5 ) 
real  straps ( 8 ) ,sumcnt,srpm,sload,spt(5) 
common/ data/tmps, gas, load, rpm , pt 

common/gas 1/pgdif f , d2 , b2 ,pgabs , wh2 ,gcf m, hw2 ,densg , 
4     viscg ,rd2  ,cd2 

common/gas 2/ ft em p,ch 4 ,co2 

common/air/ pad if f ,paabs ,d,b,airtmp,wh,acfm,hw,densa, 
&     visca,rd,cd 


read  input  file  from  device  1 ;  write  to  output  device  7 
character*40  fin,fout 
wrlte(6 ,• )iargc( ) 

if  (iargc()  .It.  2)  stop  'usage  engine  input  output' 
call  getarg( 1 ,fin) 
call  getarg ( 2  ,  f out  ) 
open( 1 ,file=fin) 
open(7,file=fout ) 
pi  =  3.141592654 

ratio  of  air  nozzle  diam.  to  plenum  intake  ^infinity 
b  =  0.0 

ratio  of  fuel  nozzle  diam.  to  fuel  tube  diam. 
b2=0.25 

air  nozzle  diam.  (in) 
d=1 .25 

fuel  nozzle  diam.  (in) 
d2r0.5 


enter  data  from  terminal 


write ( 6 ,*)'  enter  the  percent  methane  (O.xx)' 

read ( 5  ,• )frao 

write(6  ,•) 'enter  the  initial  voltages  on 
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press,  trans.  1  2  3  4  5  * 
read(5,f )  pterr! tpterr2 ,pterr3 . pterr4 ,pterr5 


call  fsetre(1) 


o   initialize  analog  channels 
1     sumcnt=0. 

do  10  i«1,8 
stmps ( i ) =0 

10  oontinue 
srpm=0 
aload=0 

do  11  1=1,5 
spt(i)=0 

11  oontinue 


100     oall  getdat 

if  (eof  ,eq.  0  .and.  count 
&        pultim  .eq.  15)  then 

srpm  =  srpm-t-rpm 
sload=sload+load 
do  12  i  =  1  ,8 
3tmps(i)=stmps(i)+tmps(i) 

12  continue 

do  13  i=1,5 
spt(i)=spt(i)+pt(i) 

13  oontinue 
8umcnt  =  sumcnt  + 1  . 
endif 


eq.  15  .and 


if  (last  .eq.  0)  go  to  100 

if  (sumcnt  .ge.  2)  then 
srpm=srpm/sumcnt 
si oad=sload/ sumcnt 
do  14  i=1  ,8 
stmps(i)=stmps(i) /sumcnt 

14  continue 

do  15  i  =  1  ,5 

spt (i) =spt(i) /sumcnt 

15  continue 


c 
c 


P8l 

padiff=.04§abs(spt(1)-pterr1) 

psi 

paabs=10.§(spt(2)-pterr2)+l4.7 
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c 
c 
c 

c 

0 


psi 

paatmspaaba 

psl 

pgdiffs.04»aba((spt(5)-pterr5)) 
psl 

Pgabs=10.»(spt(4)-pterr4)+l4.7 
ch4=frac 
oo2=1 ,-frac 
degree  R 

ftemp=(stmpa(5)+273.)-9./5. 
degree  C 

airtmp=s tmps ( 4) 


call  gasflo 
oall  alrflo 


c 
c 
c 


dyno  load  (kg) 
loadm=sload§. 45359 
dyno  torque  (ft-lb) 
torque=sload§11./12 
dyno  torque  (N-m) 
torqm=torquei1 .3825 


c 
c 

c 

c 


brake  horsepower 
bhp=2.tplttorqueisrpm/33000. 
brake  kW 
bkw=bhp».7457 

brake  mean  effective  pressure  (psi) 

bmep=H821.tbhp/srpm 

brake  mean  effective  pressure  (Pa.) 

bmepm=bmep*6894.4 


pgas=pgdiff 
gasp=(paatm+pgas)* 13.6 

fuel  gas  mass(lb/hr) 

f frgas=wh2 

mass  of  methane ( lb/hr ) 

f frch4=ffrgas»frac 

mass  of  fuel  gas(kg/hr) 

ffrgsm=ffrgas»,453  59 

mass  of  me  thane ( kg/hr ) 

ffrohm=f frgsm§f rac 
mass  of  air  (lbs/hr) 

airmas=wh 
mass  of  air  (kg/hr) 

airm=wh«. 45359 
fuel  to  air  ratio  (mass  basis) 

f toa=wh2/ wh 
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c   methane  to  air  ratio 

ch4toa=ffrch4/wh 
c   methane  to  air  ratio  (volume  basis) 

ftoav=gofm»frac/acfm 


c   calculate  specific  fuel  consumption  of  gas 
c   and  methane  ( lbs-gas/bhp-hr ) 

sfcgas=ffrgas/bhp 
c      lbs-methane/bhp-hr 

sfcoh4=sfcgas*frac 
c       kg-gas/kW-hr 

sfcgsm=f frgsm/bkw 
c       kg-methane/kW-hr 

sfcohmssfcgsm^frac 

c       calculate  thermal  efficiency 

thef  f  =  bhp»  2545./ (ffrch»»«2 1297.)*  100. 

c   calculate  correction  faotor 
cf=1 . 18*29.3 139/ paatm* 
&        ((stmps(8)+273.)/298.)".5-.l8 


o   correct  power  and  specific  fuel  consumption 
c   for  values  on  the  torque  curve 

cf bhp=bhp*cf 
cf bkw=bkw*cf 

osfcg=sfcgas*cf*bhp/cfbhp 
csfcch  =  sfcch4*cf,,bhp/cfbhp 
csfcgm=sfcgsmfcf •bkw/cf bkw 
csfcom=sf cohm^of •bkw/cf bkw 


write(7t500)yrfmo,day,hr,min,sec 
write(7i501)stmps( 1) ,stmps(2) ,stmps(3) , stmps( 4)  , 

&  stmps ( 5) , stmps( 6 )  ,stmps ( 8) 

write(7,502)loadm,sload,srpm,torqm, torque, bkw , 
&bhp,bmepm,bmep,ffrgsm,ffrgas,ffrohm,ffrchi|,airmf 
&ch4toa, f toav , sf cgsm , sf cgas , sf cohm, sf oohM , 
4thef f ,cf ,of bkw,cf bhp , csf ogm, osfcg,csfcom,osfcch 


c 
c 


500  format(6(i2,1x)) 

501  format( ' inlet  water  temperature ', 1  Ox , f 6 .2  ,/ , 


n 


4'outlet  wate 
&  oil  tem 
&  fl8x,f6 
&  'exhaus 
&      air  tem 

502  format( • load 
&  'rpm'  ,2 
2x  , ' (  '  , 
f6.2,2x 
pressur 
'gas  fu 
•<»,f6. 
(lb/hr) 
•air  ma 
ratio  m 
f4.2,') 
(lb/bhp 
'ch4  sp 
(lb/bhp 
•  therma 
•correct 
/  , '  corr 
consumpt 

f6.4f')« 

kg/bkw- 
,1x,f6.4 
endif 
if  (eof  .eq. 


r  temperatu 
perature ' , 1 
.2,/, 'fuel 
t  gas  tempe 
perature1 , 1 

-kg  (lb)', 2 
9x,f6.2,/,* 
f6. 2, •)••/• 

,'C,f6.2,» 
e  N/m2  (psi 
el  flow  rat 
*,«)•, /,»oh 

'  ,24x,f6.4, 
88  per  hour 
ass  basis  ( 
'  ,/  ,  'gas  sp 
-hr) ' ,5x,f6 
ecific  fuel 
-hr) ' ,5x,f6 
1  efficienc 
ion  factor 
ected  brake 
ion  kg/bkw- 
,/  , ' correct 
hr  (lb/bhp- 
,2x,'( »,f6. 

0)  go  to  1 


re» ,9x,f6 .2/ , 'engine 
1x,f6 .2 ,/ , 'air  temperature' 
temperature' , 17* , f 6 .2 ,/ » 
rature'  ,  9x , f 7 . 2 , /  ,  • ambient 
0x,f6.2,/) 


1xff6.2,2x, '( • ,f6. 2, »)',/, 
torque-nm  ( f t-lb) ' , 1 6x , f 6 . 2 , 
'brake  power-bkv  (bhp)'12x, 
)',/,  'brake  mean  effective 
)  ',12x,f6.2,l>x,'(',fl».2t')',/, 
e  kg/hr  ( lb/hr ) • , 24x , f 6 . 4 ,2x , 
4  fuel  flow  rate  kg/hr 
2x,'( «,f6. 4, ')',/, 

kg/hr' ,5x,f7. 4,/, 'fuel  to  air 
volume  basis)  •  ,5x,f4.2,5x, • ( • , 
ecific  fuel  consumption  kg/bkw-hr 
.4,2x,'(',f6.4,')\/, 

consumption  kg/bkw-hr 
.4,2x,'(',f6. 4, ')',/, 
y',15x,f6.2f'J',/f 
for  torque  curve  values *,  1 3x , f6 .2 , 

power  bkw  ( bhp )  '  ,2x , f 6 . 2 , 
hr  (lb/bhp-hr)  '  ,1x,f6.4,2x, •( • , 
ed  ch4  specific  fuel  consumption 
hr)  • 
4, ')',///) 


stop 
end 


subroutine  getdat 

inolude  'flog.h' 

real  tmps ( 8 ) , gas , load , rpm , pt ( 5 ) 

common/data/tmps ,gas, load, rpm, pt 


200   call  fgetre 
if  (eof  .ne. 


0)  return 


calculate  copper-constantan  temperatures 
do  210  i=1 ,5 

analog(i) =analog(i)-analog( 16 ) 
analog(8)=analog(8)-analog( 16) 
vref rfttov(analog(8)»20. ) 
do  220  i=1 ,5 

tmps(i)rfvtot(analog(i)+vref) 
tmps(8)=fvtot(vref) 


210 


220 


M8 


240 


calculate  iron-constantan  temperature 
analog(6)=(analog(6)-analog( 16) )• 1000. 
tmps(6)=tmps(8)+((analog(6)»32.72)-32.) 

•5./9. 
rpm=(analog(9)-analog(7) )• 1000. 
load=(analog( 10)-analog(7))f100. 
gas  =  pulse(  1 ) 
do  240  1=1 ,5 

pt(i)ranalog(i+10)-analog(7) 
continue 
return 
end 


c   subroutine  gasflo 


c 
c 


calculates  gas  fuel  flow  from  the  differential 
pressure  across  an  asme  nozzle. 


c  argument(s)  required  from  the  calling  routine 

c  (located  in  common  statement) 
c      pgdiff 
c      d2 
c      b2 
c      pgabs 


c   argument(s)  supplied  to  the  calling  routine 


c 

(located  in  common 

statement ) 

c 

wh2 

c 

gcfm 

c 

hw2 

c 

densg 

c 
c 

viscg 
rd2 

c 

cd2 

calculations 

subroutine  gasflo 

include  'flog.h1 

real  kf , kchl ,koo2 

real  lhvch4 , lhvco2 

common/gas 1/pgd if f , d2 , b2 , pgabs , wh2 ,gcfm,hw2 , 

densg, viscg , rd2  ,  cd2 
common/ga82/ftemp,cht,co2 


c    nozzle  discharge  coeff.  initially  assumed  to  be  0.99 

od2  =  .99 
c    nozzle  thermal  expansion  factor  is  assumed  to  be  1.0 
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c 
c 
c 
c 
c 


c 
c 

c 
c 


f  =  1  .0 

velocity  of  approach  factor  is  assumed  to  be  1.0 

f1  =  1.0 

d2  is  nozzle  diameter  in  inches. 

kf  is  gas  specific  heat  ratio  cp/cv  of  fuel 

rf2  is  fuel  gas  constant  ( f t-lbf /lbm-deg  r) 

b2  is  ratio  of  nozzle  dia.  to  intake  tube  dia 

tfuelr  =  ftemp 
pi  =  3.141592654 

calculate  the  fuel  molecular  weight 

tmwt  =  ch4«l6.  +  co2»i»4. 

partial  molecular  wt  of  fuel  due  to  oh4 

zch4=ch4»l6./tmwt 

partial  molecular  wt  of  fuel  due  to  co2 

zco2=co2#44./tmwt 

gas  constant  for  ch4  ( f t-lbf /lbm-deg  r) 

rchJ»  =  96.33 

gas  constant  for  co2  ( f t-lbf /lbm-deg  r) 

roo2=35.11 


c   calculate  fuel  gas  constant 
rf2=rch4»zoh4+roo2*zco2 

c      specific  heat  ratio  (op/cv)  of  ch4 

kch4s1.31 
c       specific  heat  ratio  (cp/cv)  of  oo2 
kco2=1 .31 

c   calculate  specific  heat  ratio  of  the  fuel 

kf rkch4§zoh4+rco2#zoo2 

c      viscosity  of  ohH  (lbf/hr-ft  at  77f) 

vch4s.0259 
c      viscosity  of  oo2  (lbf/hr-ft  at  77f) 

vco2=.0348 
c   calculate  fuel  viscosity( lbf/hr-ft  at  77f) 

viscgrvch4*zoh4+vco2*zco2 

c   convert  viscosity  to  (lbm/ft-s) 

c   multiply  by  32.2( lbm-f t/lbf-s-s)  &  divide 

c   by  32.2(ft/s-s)4  divide  by  3600(s/hr) 

c         (lbm/ft-o) 

visog=visog/3600. 
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c 
c 


lower  heating  value  for  ch4  (btu/lbm) 

lhvch4=21500. 

lower  heating  value  for  oo2  (btu/lbm) 

lhvco2=0 . 


c   calculate  lower  heating  value  of  fuel  (btu/lbm) 
lhv  =  zch4*lhvch4  +  zco2,,lhvco2 


c 
c 
c 


& 

c( 

c 

c 


calculate  lower  heating  value  of  fuel  at 
std.  temp,  of  77  deg  f  4  std.  press,  of 
14.7  psi  (btu/cu.ft.) 

Ihvv=(pgabs«l44/tfuelr)«(lhvch4«zch4/rchi» 
+lhvco2»zco2/rco2)«(537./tfuelr) 


calculation  of  gas  flow 


pressure  differential  across  nozzle( in.h2o ) 
hw2  =  pgdiff»1728.0/62.4 
pressure  ratio  at  nozzle 
rf1  =  ( pgabs-pgdif f )/pgabs 


c   expansion  factor  for  nozzle 

yr(kf/(kf-l .)»(1 .-rf 1«»((kf-l . )/kf ))/ 
&      (1.-rf1))«.5 

o      density  of  fuel  (lbm/cu.ft.) 

densg  =  pgabs" l44.0/(rf 2ftfuelr ) 

c   mass  flow  rate  of  fuel  (lbs/hr)  aoross  nozzle 
c   calculated  by  iteration  starting  with  an  assumed 
c   coefficient  of  discharge 
100    wh2  =  359.0«cd2»f»(d2»d2)»y«f 1« (hw2»densg )•• . 5 


c 
c 


c 
c 


convert  mass  flow  rate  of  fuel  to  linear  flow  rate 
(ft/s) 

vbar  =  ((wh2/3600.0)/densg)»(576.0/(pi»d2«d2)) 

reynold's  number 

rd2  =  ( vbar* (d2/ 1 2 .0) *densg ) /viscg 


o   recalculate  coefficient  of  discharge  from 
c   reynold*s  number 

ccd  =  .19436+( .152884»(alog(rd2)))-(9.7785e-03 
4       •(alog(rd2))»«2)+(2.0  903e-04»(alog(rd2))»«3) 

o   compare  calculated  cd2  with  assumed  cd2 
x1  s  abs(od2-ccd) 
if(x1  .le.  1.0e-4)  go  to  200 
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cd2  =  cod 

go  to  100 
200    cd2  =  cod 

wh2  =  359.0»cd2»f»(d2»d2)»y»f1» 
4      (hw2»densg)»».5 
o      gas  flow  rate(cfm) 

gofm  =  wh2/(60.0tdensg) 

return 

end 


c   subroutine  airflo 


c  looated  in  the  common  statement 

c  (d  b  paabs  padiff  airtmp) 

c  argument(s)  supplied  to  the  calling  routine 

c       wh 

c       acfm 

c       hw 

c      densa 

c      visca 

c      rd 

c      cd 


subroutine  airflo 
include  'flog.h' 

common/air/ padiff, paabs ,  d ,  b ,  airtmp , wh, 
4      acfm,  hw  ,  densa  ,  visca  ,  rd  ,  cd 
real  k 

c  nozzle  discharge  coeff.  initially  assumed  to  be  0.99 

od  =  .99 
c   nozzle  thermal  expansion  factor  is  assumed  to  be  1 .0 

f  =  1  .0 
c      velooity  of  approach  factor  is  assumed  to  be  1.0 

f1  =  1.0 
c       d  is  nozzle  diameter  in  inches, 
c       k  is  air  specific  heat  ratio  op/cv 

k  =  1  .4 
c       r2  is  gas  constant  ( f t-lbf /lbm-deg  r) 

r2  =  53.31 
c     b  is  ratio  of  nozzle  diameter  to  intake  tube  diameter 
o       if  intake  diameter  equals  infinity,  then  b=0.0 
o       ambient  temp,  in  deg.  rankine 

taambr  =  ( airtmp* ( 9.0/5 .0 ) +32.0) +459.69 

pi  =  3.141592654 

paatmspaabs 


calculation  of  air  flow 
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c   abs.  viscosity  of  air  (lbm/ft-s)  from  air  temp. 
c   between  450r  and  576r 

visoa  =  (-4.991211+(2.852927§taambr) 
&       •(  .  101l469e-0  2*taaabr**2))/1 .Oe+8 

c       pressure  differential  across  nozzle( in. h2o ) 
hw  =  padiff«1728.0/62.4 

c       pressure  ratio  at  nozzle 
r1  =  ( paatm-padif f )/paatm 

c   expansion  factor  for  nozzle 

y=(k/(k-1.)«(1.-r1»«((k-1.)/k))/(1.-r1))««.5 

o      density  of  air  (lbm/cu.ft.) 

densa  =  paatm* 1 44.0/(r2*taambr ) 

c   mass  flow  rate  of  air  (lbs/hr)  across  nozzle 

c   calculated  by  iteration  starting  with  an  assumed 

c   coefficient  of  disoharge 

100    wh  =  359.0»cd«f»(d«2)«y«f 1» (hw»densa )•• . 5 

c   convert  mass  flow  rate  of  air  to 
c    linear  flow  rate  (ft/s) 

vbar  =  ((wh/3600.0)/densa)«(576.0/(pi«d*d)) 

c       reynold's  number 

rd  =  ( vbar*(d/12.0)§densa)/visca 

c   recalculate  coefficient  of  discharge 
c   from  reynold's  number 

ccd  =  .19436  +  (  .15288i»»(alog(rd)))-(9.7785e-03 
&      •(alog(rd))"2)  +  (2.0903e-04»(alog(rd))»«3) 

c   compare  calculated  cd  with  assumed  cd 

x1  =  abs(od-ccd) 

if(x1  .le.  1.0e-4)  go  to  200 

cd  =  cod 

go  to  100 
200    cd  =  ccd 

wh  =  359.0»cd»f«(d»d)»y»f 1« (hw»densa)»»  .  5 
o       density  of  air ( lb/cu . f t . ) 

densa  =  paatm* 1 44 .0/ (r2*taambr ) 
c       air  flow  rate(cfm) 

acfm  ■  wh/( 60 .0*densa) 

return 

end 
c 
c   FLOG.H 
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c   Header  file  for  sym  data  logger  files 

c 

c   Michael  D.  Sohwarz 

c   May  5,  1982 

c 

character^  file 

integer  unit,  errflg,  eof,  lineno,  havlin 

integer  id,  switch,  nalog,  ndig,  npul,  tsmp, 

integer  tavg,  tsav 

integer  nyr,  nmo,  nday,  nhr,  nmin,  nsec,  noount 


integer  tavg,  tsav 

Lnteger  nyr,  nmo,  nday,  nhr,  nmin,  nsec,  ncour 

;ommon  /head/  file,  unit,  errflg,  eof,  lineno 
:ommon  /head/  havlin, id,  switch,  nalog,  ndig 
common  /head/  npul,  tsmp,  tavg,  tsav, nyr,  nmo 
jommon  /head/  nday,  nhr,  nmin,  nsec,  noount 

Lnteger  first,  last,  yr ,  mo,  day,  hr,  min,  sec 
Lnteger  count 


integer  first,  last,  yr ,  mo,  day,  hr,  min,  sec 

integer  count 

real  analog(32),  pulse(2) 

integer  dig(8),  pultim 

common  /rec/  first,  last,  yr,  mo,  day,  hr ,  min 
common  /rec/  sec,  count , analog ,  dig,  pultim,  pulse 

FGETREC.F 

getrec   get  a  record  from  the  file  and  place  in  the 

rec  sturct. 


Michael  Sohwarz 
c   May  5,  1982 
c 


subroutine  f setre( nunit ) 
inolude  'flog.h' 

integer  nunit 
lineno  =  0 
havlin  =  0 
unit  =  nunit 
eof  =  0 
errflg  =  0 
end 

subroutine  fgetre 
inolude  •flog.h' 

Integer  i 

character#1  errchr 
character*80  linein 

if  (eof  .ne.  0)  then 
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stop  'Getrec  reentered  after  EOF' 
endlf 


=  1 


yr  =  nyr 
mo  =  nmo 
day  =  nday 
hr  s  nhr 
min  =  nmin 
sec  =  nsec 
count  r  ncount 


if  (nalog  .ne.  0)  then 
do  11  j=1 , nalog, 6 
read(unit,  550)  linein 
do  12  k=8,80,12 

if  (linein(lc:k).eq.,E«)  linein(k:k)= »e« 
12     continue 

read  (linein,  520,  end=90,  errr80) 
&    (analog(i),  i  =  J ,min0( nalog  ,  J  +  5)) 
11     oontinue 

(6,520)  (analog(i),  i=1, nalog) 

lineno  =  lineno  +  (nalog  +  5)/6 
endif 

if  (ndig  .ne.  0)  then 

read  (unit,  530,  end=90,  err=80)  (dig(i),  i=i,ndig) 
(6,530)  (dig(i),  i=1,ndig) 

lineno  =  lineno  +  (ndig  +  9)/10 

endif 

if  (npul  .ne.  0)  then 


55 


read  (unit,  510,  end=90,  err=80)  pultim, 
&    (pulsed),  l  =  1,npul) 
o   (6,540)  pultim,  (pulse(i),  l=1,npul) 

llneno  =  llneno  +  1 

endif 

llneno  =  llneno  +  1 

read  (unit,  510,  end=90,  err=85)  nyr,  errchr,  nmo , 
&   nday,  nhr,  nmin,  nsec,  noount 
c   (6,510)  nyr,  errchr,  nmo, 
o&       nday,  nhr,  nmin,  nsec,  noount 

if  (errchr. ne. '/' .and . errchr. ne. '  ')  goto  85 

havlin  =  1 

eof  r  0 

return 

o    read  error  handling  code, 
c  unexpected  error 

80     write  (6,600)  lineno 
havlin  =  0 
goto  1 

c  this  occurs  normally  at  the  end  of  a  cassette  file 

85     last  =  1 

havlin  =  0 
c   (6,620)  lineno 
return 

o    read  eof  handing  code, 
c  unexpected  end  of  file 

90     write  (6,610)  lineno 

c  normal  end  of  file 

99     eof  =  1 
last  =  1 
c   (6,630)  lineno 
return 

c   FORMAT  STATEMENTS 

500    format(a7,8i4) 

510  format(i2,a1 ,5(i2,1x)  ,15) 

520  format(6(e11 .4,1x)  ) 

530  format(8(i5,1x)  ) 

540  format(i5,2f 10.0) 

550  format(a80) 

600    format( 'GETREO  unexpected  error    lineno:', 14) 
610    format(  'GETREO  unexpected  end  of  file 
&        lineno: ',14) 
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c 
& 
c 

& 

c 
c 
c 


620  format( 'GETREO  Hit  end  in  data  file 

lineno: • ,i4) 
630  format( 'GETREO  Hit  end  of  file  in  data  file 
lineno: » fi4) 
end 


FTHERMO.F 

function  fmtot(mv) 
real  iv,  fmtot 

fmtot  =  ((  (-0.77)  •  mv)  +  26.0  )  •  mv 

return 

end 


function  fttom(t) 
real  t,  fttom,  sqrt 

c   /•  mv  =  (26  -  sqrt(676  -  3.08  •  t))/1.54   »/ 
fttom  =  16.88311688  -  sqrt ( 285 .0396356889 

4       -  1.2987012987  •  t) 
return 
end 

function  fvtot(v) 
real  v,  fvtot 

fvtot  =  ((  (-0.77e+6)  •  v)  ♦  26.0e+3  )  •  v 
end 


c 
& 


function  fttov(t) 
real  t,  fttov,  sqrt 

/•  mv  =  (26  -  sqrt(676  -  3.08  •  t )  )/1  .54   •/ 

fttov  =  16.88311688C-3  -  sqrt ( 285 .0396357e-6 

-  1 .298701299e-6  ■  t) 

return 

end 


c   FTRANS.F 

subroutine  ftrans 
include  •flog.h1 


10 


do  10  i= 1 ,nalog 

if  (i  .ne.  7)  then 

analog(i)  r  analog(i)  -  analog(7) 

endif 

continue 


analog(8)  =  analog(8)  *  20 

oompen  ■  f t tov ( analog ( 8  )  ) 
write(6,f)  compen,  analog 
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do  20  i  =  1  ,nalog 


20 


if(i.ne.8 

analog(i) 

endif 

oontinue 

return 

end 


and . i . ne.7 .and . i . ne. 1 1 .and . i 
=  fvtot (analog(i)  +  compen) 


ne. 16 )then 


c 
c 
c 
c 


LOGVAR.F 
calculate 


the  variance  of  the  analog  data  channels 


lnolude  'flog.h' 

real  sum(32),  sumsq(32) 
character§80  fname 

call  getarg(1,  fname,  80) 
vrite(6,f)  fname 
call  setfil(1,  fname) 
call  fsetre  (1) 


100 


101 


10 


600 


do  101  1=1,32 
sum(l)  =  0 
sumsq(l)  =  0 
continue 
n  =  0 

call  fgetre 

if  (eof.ne.0)  then 

stop  'end  of  data' 

endif 

oall  ftrans 

do  10  i=1 ,nalog 

sum(i)  =  sum(i)  ♦  analog(l) 

aumsq(i)  =  sumsq(i)  ♦  analog ( I ) *anal og ( i ) 

continue 

n  =  n+1 

if  (last.eq.0)  goto  1 

if  ( n.eq . 1 )  then 

div  =  1 

else 

div  =  n-1 

endif 
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20 


xaq  =  -xsq 

endif 

sumsq(i)  =  sqrt (xsq/div ) 

if  ( sumsq ( i  )  .  ne . 0 . 0 )  then 

varian  =  1 00  .0§sumsq ( i ) /sum( i) 

else 

varian  =  0 

endif 

vrite(6,«) 

continue 

write(6,») 

goto  100 

end 


i,  sum(i),  sumsq(i),  varian 


c 

c 

c   TEST.F 

#inolude  "flog 


h" 


call  setrec(1) 

call  getrec ( ) 

write(6,§)  file,  unit,  errflg,  eof,  lineno 
write(6,§)  havlin,id,  switoh,  nalog,  ndig,  npul 
write(6,f)  tamp,  tavg,  tsav,nyr,  nmo ,  nday,  nhr 
vrite(6,f)  nmin,  nsec,  ncount 


write(6,§)  first,  last,  yr 


mo 


day,  br,  min,  sec 


write(6,f)  count , analog ,  dig,  pultim,  pulse 

if  (eof  .eq.  0)  goto  1 
vrite(6,«)  'got  EOF' 
call  getrec() 
vrite(6,«)  •oall  2' 
call  getrec ( ) 


stop  'from  main' 
end 
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ABSTRACT 

Anaerobic  digester  gas  from  two  facilities  were  individually 
evaluated  for  performance  in  an  internal  combustion  engine.  The 
first  test  series  used  digester  gas  produced  from  beef  manure 
while  the  second  series  of  tests  measured  the  performance  on  a 
swine  derived  gas. 

The  beef  based  gas,  produced  at  a  research  digester,  had  an 
average  methane  content  of  approximately  53  percent  although  week 
to  week  variations  in  the  methane  content  were  sometimes  as  great 
as  ten  percent.  The  swine  based  gas,  from  a  more  steady  state 
commercial  operation,  had  a  6*1  percent  methane  content  with  vari- 
ation of  less  than  two  percent. 

Gas  chromotography  tests  on  the  swine  based  gas  revealed 
that  in  addition  to  the  primary  methane  and  carbon  dioxide  gas 
constituents,  that  approximately  one-half  percent  of  the  gas  was 
hydrogen  sulfide  which  may  be  a  cause  of  premature  engine 
failure.  The  gas  chromotography  also  revealed  that  the  filters 
used  on  the  system  to  reduce  the  hydrogen  sulfide  concentration 
made  insignificant,  if  any,  reduction  of  the  hydrogen  sulfide 
content . 

Engine  tests  were  conducted  using  a  four  cylinder,  spark 
ignition  engine  with  an  8  to  1  compression  ratio.  The  200-  hour 
engine  test  cycle  was  the  alternative  fuel  screening  test  recom- 
mended by  the  Engine  Manufacturer's  Association.  Engine  perfor- 
mance,  lubricating   oil   deterioration,   and   engine   wear    and 


deposits  were  measured  or  monitored  to  determine  problems  encoun- 
tered with  the  use  of  the  fuel. 

Engine  performance  indicated  that  the  torque  and  power  out- 
put of  the  engine  fueled  with  digester  gas  was  from  80  to  95  per- 
cent of  the  outputs  of  the  engine  fueled  with  natural  gas.  The 
output  depended  on  the  methane  concentration  in  the  fuel,  the 
amount  of  water  vapor  present,  and  the  gas  supply  pressure. 

Modifications  required  for  satisfactory  engine  operation 
included  a  higher  fuel  supply  pressure  than  that  needed  for 
natural  gas.  The  engine  performance  was  the  best  with  gas  pres- 
sures of  7  to  12  kPa  (35  to  50  inches  of  water).  An  increase  in 
the  fuel  to  air  mixture  and  increased  spark  advance  were  also 
required . 

Hear  metal  concentration  levels  in  the  lubricating  oil 
revealed  no  excessive  engine  wear.  Total  Base  Number  of  the  oil 
remained  high  with  oil  change  intervals  of  100  hours.  Main  bear- 
ing clearances  after  MOO  hours  of  operation  on  the  digester  gas 
were  well  within  the  tolerances  of  new  bearings.  Visual  inspec- 
tion of  other  engine  components  revealed  no  excessive  wear. 

A  problem  found  in  the  first  test  was  excessive  amounts  of 
water  vapor  in  the  gas.  Deposits  on  the  engine  components  were 
minimal  after  the  tests  on  the  beef  digester  gas,  but  after 
operation  at  the  swine  facility,  the  deposits  were  significant, 
engine  failure  could  have  ocourred  within  a  short  time  period. 


